P-wave Cooper pair splitting 
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Splitting of Cooper pairs has recently been realized experimentally for s-wave Cooper pairs. A 
split Cooper pair represents an entangled two-electron pair state which has possible application in 
on-chip quantum computation. Likewise the spin-activity of interfaces in nanoscale tunnel junctions 
has been investigated theoretically and experimentally in recent years. However, the possible impli- 
cations of spin-active interfaces in Cooper pair splitters so far have not been investigated. 
We analyse the current and the cross correlation of currents in a superconductor ferromagnet beam- 
splitter including spin-active scattering. Using the Hamiltonian formalism we calculate the cumulant 
generating function of charge transfer. As a first step, we discuss characteristics of the conductance 
for crossed Andreev reflection in superconductor ferromagnet beamsplitters with s-wave and p-wave 
superconductors and no spin-active scattering. In a second step, we consider spin-active scattering 
and show how to realize p-wave splitting only using a s-wave superconductor via the process of spin- 
flipped crossed Andreev reflection. We present results for the conductance and cross correlations. 
Spin-activity of interfaces in Cooper pair splitters allows for new features in ordinary s-wave Cooper 
pair splitters, that can otherwise only be realised by using p-wave superconductors. In particular it 
provides access to Bell states different from the typical spin singlet state. 

PACS numbers: 



I. INTRODUCTION 



Solid state entanglers represent the electronic ana- 
logues to the optical setups used for Bell inequality 
tests. Most setups propose a superconductor as a source 
of spin-entangled s-wave Cooper pairs 1 - - — . These are 
transferred to spatially separated leads by the process of 
Crossed Andreev Reflection (CAR) 4 .^. Different setups 
have been considered in order to achieve CAR without 
being dominated by elastic cotunnelling or ordinary 
Andreev Reflection (AR) processes^—. S-wave Cooper 
pairs are entangled in energy and spin space. Therefore, 
one may either filter the electrons of a Cooper pair 
in spin space (using ferromagnets 1 ^— or Luttinger 
liquids^), or in energy space (using quantum dots 1 ^, 
coupling to an electromagnetic modeM, an appropriate 
voltage biasi^— or ac-biasi^). Filtering using quantum 
dots with large onsite interaction has been successfully 
realized in experiment^— and a non-local resistance has 
been measured as characteristics of Cooper pair splitting. 
Moreover, a positive cross correlation measured in a 
superconductor and normal-metal three-terminal device 
gave compelling evidence for CAR 23 . However, s-wave 
Cooper pairs only give access to one of the Bell states, 
namely = 1/V2 (| f)i ® U>2 - U)i ® | t>a), where 
the subscripts 1 and 2 refer to two normal leads. P-wave 
Cooper pairs give access to the other Bell states, espe- 
cially those like = l/y/2 (U)i<g> \±) 2 ± | t)i ® 1 1) 2 ), 
involving fully spin-polarized combinations of the two 
electrons. Thus, a p-wave Cooper pair splitter represents 
the essential counterpart to s-wave Cooper pair splitters 
as on-chip sources of spin-entangled Einstein-Podolsky- 
Rosen (EPR) electron pairs^— . 

In this paper we approach the problem of p-wave Cooper 



pair splitting in two steps. First, we show that p-wave 
splitting may easily be identified in a hybrid junction 
between a superconductor and two ferromagmets (setup 
shown in Fig. [l}. However, p-wave superconductors that 
can be easily handled in quantum transport experiments 
are presently not available. Therefore, in a second step 
we show how p-wave splitting can be realized without 
using a p-wave superconductor. Indeed, the previous 
works on the charge transfer statistics of superconductor 
ferromagnet beamsplitter s] 10 ' 27 ' 28 neglected the recently 
predicted22r— and observe d 32 ' 33 effects of the interface 
on charge transfer. By including a better description 
of the interface one can describe the shared triplet 
pairs formation between the ferromagnets. We calculate 
the Full Counting Statistics (FCS) which allows us to 
identify the corresponding charge transfer process and 
to calculate the cross correlation as an experimentally 
observable quantity^. 



II. RESULTS AND DISCUSSION 



A. Superconductor Ferromagnet Beamsplitters 



Splitting of spin-polarized p-wave Cooper pairs can 
easily be identified in the conductance. From the result 
in2I we find the generalization of Beenakker's formula 34 
for zero-bias conductance of a beamsplitter realized by a 
resonant level between a superconductor and two ferro- 
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magnets 

GcA R (cri,a 2 ) = [4r 2 r2(l + (Tl P 1 )(l + a 2 P 2 )] 

x [A3 + r| + r 2 (2 + o-iPi + <j 2 p 2 ) 2 } - 1 

x [A 2 + r| + T 2 (2 - o-iPx - a 2 P 2 ) 2 ] "(1) 

where T is the tunnel rate through the two barriers be- 
tween the quantum dot and the ferromagnet (both are 
assumed to be equal) and V g is the tunnel rate between 
the superconductor and the quantum dot. A^ refers to 
the energy of the resonant level. Pi and P 2 are the (par- 
allel) polarisations of the ferromagnets and <j\ and <r 2 are 
the spins of the electrons in a Cooper pair. 
In usual s-wave superconductivity the spin directions 
obey o\ = — er 2 , and we may maximize CAR by choos- 
ing Pi = 1 = — P 2 (or vice versa)**. However, if we 
choose the polarisations of the ferromagnets to be equal 
(i.e. Pi = 1 = P 2 ) Gcar in Equation ([]} becomes zero 
(see Fig. [IJ. The situation is reversed if we introduce 
a spin-polarized p-wave superconductor so that a\ = cr 2 . 
Now splitting is maximized if Pi = P 2 . Of course now for 
antiparallel polarisation the current is blocked (see Fig. 
[T]). Therefore p-wave splitting is easily identified in the 
crossed conductance. The assumption of maximal polar- 
isation is not easy to realize in experiment 3 ^ but we use 
it to illustrate the argument. 

In spite of the apparent simplicity of this consideration 
a serious problem remains: p-wave superconductors are 
very rare and materials such as Sr 2 Ru04 are hard to 
handle^. 

However, recent theoretical and experimental progress 
showed that the treatment of superconductor - ferromag- 
net interfaces requires special care with respect to the 
exact form of the interface and the associated spin-active 
nature of tunnelling 3 -^ In the rest of this work we want 
to show how a spin-active interface can be used in exactly 
the same way as a p-wave superconductor and therefore 
allows to generate the other Bell states. In order to ac- 
complish this goal we will show how the transport charac- 
teristics of p-wave splitting can be imitated by spin-active 
scattering. 

B. Cumulant generating function with spin-active 
scattering 

To identify the separate charge transfer processes and 
to evaluate the conductances and cross correlations we 
calculate the FCS of charge transfe r 38 ' 39 using the gen- 
eralized Keldysh technique^2r— . The aim is to calcu- 
late the Cumulant Generating Function (CGF) x(A) of 
the probability distribution P(q) to transfer q units of 
charge during a (long) measurement time r: x(A) = 
E 9 e iXq P{q). Partial derivat ives of this function provide 
direct access to the cumulants (irreducible moments) of 
P{q). 

We model the superconductor ferromagnet beamsplitter 



by two ferromagnets PI and P2 which are tunnel cou- 
pled to a resonant level, which is the simplest model of a 
quantum dot. The superconductor is also coupled to the 
resonant level via Hti and Ht2 which takes into account 
the interface effects. The result is the Hamiltonian 

H = Hpi + Hp2 + Htfi + Htfi + Hs (2) 
+P T i + P T2 + H d , where (3) 

H d = ^A d J+d CT , (4) 

a 

represents the resonant level. Throughout the rest of 
exposition we use units in which e = h = ks = 1 and 
restrict ourselves to a quantum dot on resonance A^ = 0. 
Ferromagnetic electrodes are described in the language of 
electron field operators ^k,a.a, where a = P1,P2 using 
the Stoner model with an exchange energy h ex as in^ 3 - 

-h ex , a X)(*fc,a,t**.«.t - *fc,«a*k,<a)- ( 5 ) 
k 

Consequently they can be described as fermionic con- 
tinua with a spin-dependent DOS po,a,a = Po,a(l + cP a ), 
where P a is the polarisation. The superconductor is de- 
scribed using the ordinary s-wave BCS Hamiltonian 

fc,tr 

+ A Ys^ts^Usi + 9-kSi*kin). (6) 
k 

We express the energies of the dot and the reservoirs 
relative to the superconductor chemical potential^ so 
that [is — and V a — /is — Ha = —fJ-a is the chem- 
ical potential of the ferromagnets. The quasiparticle 
density of states in the superconductor has the form: 
Ps = PosM/Vw 2 — A 2 > where pos is constant in the 
wide-band limit. Tunnelling between the superconduc- 
tor and the quantum dot is given by^ 

ffrs = X)7s[<£*s )CT (z = 0) + H.c.] ! (7) 

a 

where 7s is the corresponding tunnelling amplitude be- 
tween the dot and the superconductor. The tunnelling is 
assumed to be local and to occur at x = in the coor- 
dinate system of the respective electrode. ^s,a{x) refers 
to the electron field operator of the superconductor in- 
troduced in Equation ([6]) in position space. 
Finally, we need a Hamiltonian approach^ for spin-active 
scattering. There are manifold effects like spin-orbit cou- 
pling, magnetic anisotropy or spin relaxation that give 
rise to spin-activity of interfaces 29 . Previous studies of 
point contacts used a scattering states description 2 ^ 3 1 1 47 
in order to introduce a spin-active scattering angle as a 
phenomenological parameter to characterise the interface 
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FIG. 1: Summary of the possible charge transfer processes in a superconductor ferromagnet beamsplitter. The superconductor 
(blue) and the two ferromagnets (red, are assumed to be fully polarized) are coupled via a quantum dot realised by an InAs 
nanowire (brown). The polarisation is indicated by green arrows. In the upper part, the situation for s-wave superconductors is 
shown. The Cooper pair may split if the two ferromagnets are antiparallely polarized. In the lower part the reversed situation 
for spin-polarized p-wave superconductors is depicted. The Cooper pair may now split only if the two ferromagnets are equally 
polarized. 





or a wave-function matching techniqu e) ' . We adopt 
the simplest possible approach and follower— by intro- 
ducing two tunnelling Hamiltonians at each of the inter- 
faces between the quantum dot and a ferromagnet 



Htfi, a = ^2%,i[<%*<*A x = °) + Hc L 
H TF 2, a = 5^7a, 2 [#¥a,-a(a: = 0)+H.a]. (8) 



Htfi. a describes normal spin-conserving tunnelling 
whereas Htf2,cx refers to the spin- hip processes 
[^ aj - a (x = 0) represents an electron in the ferromagnet 
with opposite spin] . If we take spin-active scattering into 
account this way we have 5 different tunnel couplings. In 
order to reduce the number of parameters and for clari- 
fication of the discussion of spin-active scattering effects 
we want to limit ourselves to a special constellation of 
parameters as far as spin-active scattering is concerned, 
namely 7fi,i = Jf2,i, 7fi,2 = 7F2.2- In this case we 
define operators 



\pp~\,i + 7fi,2 



which allow us to rewrite the tunnelling Hamiltonians in 
Equations and jSJ) as 

H TF , a = 5^7 o [4*a,a(a: = 0)+H.c.] ) (10) 

Hti = ^7s,i[4*s,a(z = 0)+H.c.], 

H T 2 = 537s ) a[4*s,-a(a: = 0)+H.c.]. (11) 

a 

This minimal model still reveals all the transport prop- 
erties we want to discuss here. 

In order to access the CGF, we introduce the count- 
ing fields for the leads attached to the quantum dot as 
time-dependent fictitious parameters aka counting fields 
in front of the creation and annihilation operators of the 
respective electrodes 

*sAx = 0) -> e- lXs ^'HsAx = 0), 
*a )CT (z = 0) -> e-'^W/^^^^O). (12) 

\s,i,2)(t) takes the value +(—)Xrs,i,2) on the for- 
ward/backward branch of the Keldysh contour C. The 
counting fields are only switched on during the measure- 
ment time t. In the limit of long measurement times 
r the CGF can be calculated analytically using a gen- 
eralized Green's function formalism that has previously 
been used to calculate the FCS of other quantum impu- 
rity systems^—. Following^ the CGF is given by the 
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expectation value 
1iiXsff(A, r) 



T c exp 



_rp\s 
± T1 



T, 



TF2 



T, 



12 



(13) 



where A = (Ai,A2,As) and Tc denotes the contour time 
ordering operator. T£ FV ■ ■ ■ , T T | are abbreviations for 



the tunnelling operators introduced in Equations (|10[) 
and (jlip in combination with the substitutions defined 
in Equation ([12"]) . 

For the FCS calculation in the limit of large measurement 
times (as assumed here) we follow^ - and write the above 
expression using an adiabatic potential^ which, in the 
limit of infinitely long measuring time, does not depend 
on it 



In Xsff(A) 



dtU(X) = -irU{X). 



In turn the adiabatic potential is related to the counting 
field derivatives of the T x, s introduced in Equation ([13")) 



3=S,1,2 



dU 



dX 



•3=5,1.2 




where (-)\ is defined as (-)\ :— 1/xsff(A, t)(-) with (•) 
being the ordinary expectation value with respect to the 
systems Hamiltonian in Equation (j3|) with the tunneling 
Hamiltonians rephrased using the substitution defined in 
Equation (fl~2j). The derivatives of the T A 's take the form 
of mixed and A-dependent Green's function, e.g. 



to correlation functions of the type (d^d^)\— . Conse- 
quently the full dot Green's function for the spin species 
f becomes a 4-component vector of correlation functions 
{(dfd^}x,{d t d^}\,{d^dpx,( d ^ d: t)x)- Since tne count- 
ing fields take different signs on the backward/forward 
branch of the Keldysh contour each correlation function 
is a 2x2 Keldysh matrix again. The bare dot and lead 
Green's functions are e.g. given ir* 6 ^. Since the solution 
of the Dyson equation allows to sum up all orders in the 
tunnelling couplings the result is exact and also valid at 
finite temperature as discussed also in—. 
However, the result for the average in Equation (fT"3|) is 
quite long. We want to restrict ourselves to the study of 
the relevant aspects of the CGF in view of the possbility 
of p-wave splitting only. 

Comparing our method to previous treatments of super- 
conductor hybrids we should emphasize that also in the 
quasiclassical Green's function formalism 61 calculations 
of the current exact in tunneling have been don o 62 ' 63 . 
However, there it was assumed to have two tunnel con- 
tacts to normal leads with large spatial separation. In 
our case the three leads are coupled via a quantum dot 
which corresponds to the opposite limit and involves an 
energy-dependent transmission. Therefore the setup con- 
sidered here is closer related to the studies in Refs. 58 
and 6 ^ where, however, a disordered quantum dot and 
no spin-active scattering was assumed, whereas here we 
treat the ballistic case with spin-active scattering. Disor- 
der for the case of tunnel contacts was considered using 
the Keldysh-Usadel formalism in Refs^r 6 - and 6 -^. Ad- 
ditionally we want to point out that the description of 
spin-active scattering used here and in Refsi^ 3 . and 2 ^ is 
different: in the latter works the spin-mixing angle is in- 
troduced as a phenomenological parameter whereas we 
use a second tunneling transparency to account for spin 
flips. For a comparison of both descriptions of tunnel 
contacts we refer the reader to^ 8 -. 



T, 



TFl 



dX x 



-iAi/2 



l,ff/A 



h.C. 



C. Spin-flipped Crossed Andreev Reflection 



These mixed Green's functions have to account for all 
orders of the tunneling coupling. By expanding them to 
first order in the relevant tunnel coupling they can all 
be rewritten as a product of a bare lead Green's func- 
tion (subscript 0) and an exact-in-tunneling dot Green's 
function. E.g. for the derivative considered above 

(T c d+*i l0 .) A ex 7i<#i, CT ^ CT )o<eUj}A- 

Therefore the remaining task is an exact calculation of 
the dot Green's function. Since the system's Hamilto- 
nian in Equation ([3]) is quadratic the Dyson equation 
can be solved exactly. However, as also noted in^S, cou- 
pling of the dot to the superconductor automatically 
leads to the appearance of anomalous dot Green's func- 
tions of the type (didf)x to second order in the tunnel 
coupling to the superconductor. Additionally, the spin- 
flipping tunnel contribution in Equation (1111) gives rise 



In simple superconductor- ferromagnet tunnelling junc- 
tions the presence of spin-active scattering gives rise to 
a new type of AR. In ordinary AR an electron is retrore- 
flected as a hole with opposite spin since Cooper pairs 
represent spin singlets. However, due to the spin-active 
nature of tunnelling in the setup considered here the hole 
or the electron spin can be flipped. This Spin-flipped 
Andreev Reflection (SAR) induces triplet correlations in 
the ferromagnet™ To see whether similar effects occur 
in our setup we calculate the zero-bias conductance anal- 
ogously to Equation JTJ) for the case of Pi = P2 = 1, 
meaning maximal parallel polarisation 



Gscar 



16e 2 



4r 2 r| F (4r 2 



SF 



h [r 



8r l r SF 



ier 2 F (4r 2 



r 2 

1 SF 



Ts = ^Pashsi + ls2) and r F,a = * 'Po7« = : r a g ain refer 
to the tunnel rates for ordinary single electron tunnelling. 
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= 7rpos7si7s2 describes the additional spin-flip tun- 
nelling rate at the interface. Obviously, spin-active scat- 
tering has lifted the current blocking indicated in Fig. Q] 
for a s-wave superconductor connected to two maximally 
parallelly polarised ferromagnets. Therefore this finite 
conductance is similar to the one obatined for a p-wave 
superconductor junction in Equation This indicates 
that this conductance for voltages below the gap is indeed 
associated to a Spin-flipped Crossed Andreev Reflection 
(SCAR) in which a triplet pair is transferred to the fer- 
romagnets. 

Of course one may also obtain this information from the 
CGF itself. However, the expression is complicated and 
the probability distribution of charge transfer is also not 
easy to access in an experimen t 69 ' 70 . Therefore we follow 3 
and use the cross correlation as an indication to probe 
whether the two charges of a Cooper pair transferred in 
a SCAR process are indeed transferred to the ferromag- 
nets at the same time. The presence of non-zero Gscar 
at voltages below the gap and T = in combination with 
a positive cross correlation can only be explained by a si- 
multaneous transfer of a triplet pair to the ferromagnets, 
which implies that we indeed observe p-wave splitting. 
The cross correlation can be calculated as a mixed sec- 
ond derivative of the CGF 



P 1 
1 1 ■ 



1 d 2 In xsff (A, r) 



d\\d\2 



(14) 



In Fig. [2] the result of Equation (fT4| is shown for two 
different configurations of the couplings and polarisa- 
tions. For moderate polarisation (P = 0.3) and no 
spin-active scattering we find two cases in which positive 
cross correlation can be observed in accordance with 
previous result a 27 ' 71 . First, for voltages close to the 
superconducting gap and V\ ~ — Vi CAR is strongly 
suppressed and one expects single electron transmission 
to be dominant. Nonetheless, the energy-dependent 
DOS of the superconductor leads to large transmission 
coefficients for double AR from one ferromagnet to the 
superconductor and further to the second ferromagnet. 
This process is known as Andreev reflection Enhanced 
Transmission (AET) 2 - 7 -. The second case of positive cross 
correlation can be observed for one bias voltage being 
close to zero and finite bias on the second electrode. In 
this case, CAR dominates over single electron transmis- 
sion and induces positive cross correlation^ 7 -. 

This picture changes dramatically if we go over to the 
case of full polarisation (P = 1) and finite spin-active 
scattering. Since AET relies on double AR, and con- 
sequently the spin symmetry of AR, it must dissapear 
since SAR violates spin-symmetry and it would be the 
only possible charge transfer process for a single lead for 
P = 1. However, for Vi Va a positive cross correlation 
remains. This is exactly the position where we assume 
SCAR to be dominant since V\ ps V2 means that single 
electron transfer between the ferromagnets does not 
occur. 

The effect is, of course, still observable for P < 1 but 



the polarisation should be rather strong. Spin-active 
scattering in superconductor-ferromagnet hybrids is a 
general phenomenon™ and full polarisation was just 
assumed for clarity. Therefore, we believe that SCAR 
is a generic phenomenon that should also be present in 
superconductor ferromagnet beamsplitters with tunnel 
contacts 6 -!^ 3 - or chaotic cavitie o 58 ' 64 since its origin does 
not lie in the precise form of the energy dependence of 
the transmission coefficients. 

Experiments in the direction of the above described 
proposal have already been realized 7 - 3 -. Multi-terminal 
hybrid systems with embedded quantum dots 7 - 4 - also 
based on InAs nanowire o 26 ' 75 have already been realised 
experimentally In such devices a new subgap structure 
has been observed that can be explained by SAR 68 . This 
is of special importance, since in our consideration we 
did not include effects of Coulomb interaction on the 
dot and so we should worry about a possible a suppres- 
sion of SCAR. The mean field analysis o&, however, 
revealed that SAR and therefore also SCAR should be 
observable also in presence of strong Coulomb interac- 
tion. Apart from that, we can argue that also in in- 
teracting systems characteristic resonances, as the one 
of the resonant level considered here, are present and 
have a characteristic location and width associated with 
interactions. Therefore the general scenario should be 
robust. However, one should mind that for more ex- 
tended quantum dots or nanowires disorder could play 
an important role as mentioned in££ and^ 6 -. Another 
experiment has realized a superconductor-ferromagnet- 
superconductor junction based on Al and Co electrodes 
with two closely spaced cobalt wires bridging two alu- 
minum electrodes^ 6 -. In this experiment the resistancee 
in the case of antiparallel and parallel magnetisation of 
the two wires has been measured and for low tempera- 
tures it has been found that the antiparallel arrangement 
may even have higher resistance than the parallel one 
giving reliable evidence for spin-active scattering being 
present in the device. Concerning a possible experimen- 
tal realisation using quantum dots one should mind that 
interaction should be small enough and the polarisation 
large enough to not completely suppress SCAR and the 
coupling to the leads should be of the order of the su- 
perconductor gap r, Ts, T5^«A. Nowadays, this cou- 
pling is generically obtained in experiments using InAs 
nanowires or carbon nanotubes as quantum dots 7 - 7 -. 



III. CONCLUSION 

To conclude, we have first considered supercon- 
ductor ferromagnet beamsplitters without a specific 
consideration of interface properties. We found that 
splitting of spin-polarized p-wave Cooper pairs can 
easily be identified in the conductance. However, p-wave 
superconductors that are usable in experiments are 
not available but we proposed a scheme to mimic its 
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FIG. 2: Cross correlation of currents for different parameters: The cross correlation P 12 according to Equation (|f4|) is calculated 
for two different sets of parameters. The polarisation of the ferromagnets is assumed to be equal in both cases. The left graph 
shows the result for P = 0.3, T s = 4A, F S f = 0, F F i = 0.4A, F F2 = 0.1A and T = 0.1A. The right graph is for P = 1, 
Fs = 2A, Fsf = A, Ffi = 0.05A = Ff2 and T = 0.1A. A distortion of the superconducting DOS described by a Dynes 
parameter of Fd = 0.005A has been introduced in order to circumvent numerical artifacts of the diverging superconductor 
DOS. The cross correlations are positive in the deep red regions. 



behavior by taking into account the spin-activity of 
superconductor-ferromagnet interfaces. The newly 
identified SCAR process allows to obtain split p-wave 
Cooper pairs which gives access to the Bell states other 
than |* _ ). 
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